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Pion form factors from lattice QCD with exact chiral symmetry 



1. Introduction 

A detailed study of the chiral behavior of pion observables is an important subject towards 
a deep understanding of the low-energy dynamics of QCD. The pion vector form factor Fy(q 2 ) 
has been precisely measured by experiments and the charge radius (r 2 )y can be extracted in a 
model-independent way. Although a detailed comparison of (r 2 )y between chiral perturbation 
theory (ChPT) and non-perturbative determinations in lattice QCD is very interesting, the chiral 
behavior of (r 2 )v is distorted by explicit chiral symmetry breaking with the conventional lattice 
quark actions. 

The chiral behavior of the scalar form factor Fs(q 2 ) is another interesting subject. The scalar 
radius (r 2 }s provides a determination of the low-energy constant (LEC) U, which should be com- 
pared with that from F K , and (r 2 )s has 6 times enhanced chiral logarithm compared to (r 2 )y. Since 
there are no experimental processes directly related to Fs{q 2 ), its direct determination is possible 
only through lattice QCD. It is, however, computationally very demanding to evaluate contributions 
of disconnected diagrams to Fs{q 2 ) with the conventional simulation methods. 

In this article, we present our calculation of these pion form factors in two-flavor lattice QCD. 
Chiral symmetry is exactly preserved by employing the overlap quark action ^ for a direct com- 
parison of (r)v,s with ChPT at two loops. We use the so-called all-to-all quark propagator [Q] to 
calculate Fs(q 2 ) including the disconnected contributions. We refer readers to Ref. [||] for more 
detailed description of this study. 

2. Simulation method 

We simulate QCD with two flavors of degenerate up and down quarks using the Iwasaki gauge 
action and the overlap quark action. We also introduce a topology fixing term [Q] into the lattice 
action to substantially reduce the computational cost. Our numerical simulations are carried out 
on a xN t = 16 3 x 32 lattice at a lattice spacing of a = 0. 1184(21) fm, which is fixed from the 
heavy quark potential. In the trivial topological sector Q = 0, we simulate four values of the bare up 
and down quark masses m uc i = 0.015, 0.025, 0.035 and 0.050, which cover a range of the pion mass 
290<M^[MeV] <520. Statistics are 100 independent configurations at each m U( j. We also simulate 
nontrivial topological sectors Q = — 2 and —4 to study effects of the fixed topology. Further details 
on our configuration generation are presented in Ref.[Q]. 

The matrix element (7l(p')\0r(p' — p)\k(p)) can t> e extracted from the three-point function 

C n( ? rn (At,At';p,p') = ^^£c^ rJr (x,f;Af,Ar';p,p'), (2.1) 
C K ^{x,t;At,At';p,p') = £(^(x^^ + A^ + AO^^(x^^ + A0^(x,0 t ^'- ^ ' p ' (x '- x '' ) e - ! ' p(x ''- x) ),(2.2) 

x'x" 

where &y is the vector current or scalar operator S, and G n represents an interpolating field to 
create the physical pion state. In the conventional method, one calculates the so-called point-to-all 
quark propagator Sp(x' ,x), which flows from a fixed lattice site x to any site x', by solving 

Y,D(y,x')S F (x',x) = 8 y , x , (2.3) 

x' 

where D is the Dirac operator. Then, C 7[ ^ r7t (x,t; At, A?';p,p') can be calculated by connecting the 
point-to-all propagators as shown in Fig. [j]. However, we need to solve Eq. ( p73| ) for each lattice site 
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Figure 1: Connected (left-most diagram) and disconnected (middle diagram) three point functions. The 
right-most diagram with the vacuum expectation value (VEV) (^r) contributes to Fs(Q). Each line shows 
quark propagator Sf- We use Sf(x,x!) — YsSf(x\x)^ 75 for the spectator quark. 



x to evaluate the disconnected quark loop D~ l (x,x) as well as to average C n ff Tn (x,t; At, Af';p,p') 
over the pion source location x = (x,f) as in Eq. (jO)). This needs prohibitively large CPU cost. 

This difficulty can be avoided by constructing the all-to-all quark propagator, which contains 
the quark propagating from any lattice site to any site, in an effective way. Along the strategy 
proposed in Ref. [Q], we prepare 100 low-lying modes of D for each gauge configuration, and their 
contribution to the all-to-all quark propagator is calculated exactly as 




(2.4) 



where Xk and represent k-th smallest eigenvalue and corresponding eigenvector. It is expected 
that this low-mode contribution dominates low-energy dynamics of pions. Possibly small contribu- 
tion from the remaining high-modes can be estimated stochastically by the so-called noise method 
[Q], which is not computationally intensive. We refer to Ref. [J3j] for more technical details on our 
method to calculate pion correlators using the all-to-all propagator. 

We calculate three-point functions, C % y A% and C n $ n , as well as the two-point function 

C™(A/;p) = ^££(t? n (x',t + At)t? n (x,t) f e-W-^). (2.5) 

We take 33 choices for the spatial momentum pM with |p^| < 2, which cover a region of the 
momentum transfer — 1.7 < q 2 [GeV 2 ] <0 for the three-point functions. Note that the spatial mo- 
mentum is shown in units ofln/L in this article. 



3. Determination of pion form factors 

For a precise lattice calculation of the form factors, it is advantageous to construct an appro- 
priate ratio of correlators so that some uncertainties, such as renormalization factors, cancel at least 
partially in the ratio []7|]. In this study, we calculate effective value of the vector form factor from 

F (At At' ^ - 2Mn Rv{&,At';\v\,\v'lq 2 ) 

H ' ,q) ^(IpD+^Gp'I) R v (At,At';0,0,0) ' ( j 

R (At Af'-lnl In'l n 2 \ - 1 V C^jAt ,At' 

Rv(At,At ,|p|,| P u ) - ^^^(^rtq^^pO' (3 " 2) 

where (l/Af|pMp/|)£fixed IpUp'l represents the average over momentum configurations {(p,p')} cor- 
responding to the given value of q 2 . As shown in the left panel of Fig. ^, we observe a very clear 
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Figure 2: Left panel : effective value Fv(At ,At';q 2 ) at m„</ = 0.025, which is around a quarter of the physical 
strange quark mass. Right panel : statistical fluctuation of C w y 4W (A/,A?';p,p') at (|p|, |p'|) = (y/2,0) and 
Af = At' = 7. Filled and open symbols are results with and without averaging over the source location (x,f ). 



signal of Fy (At, At';q 2 ) with the statistical accuracy of typically 3-5%. We note that the use of the 
all-to-all propagator enables us to achieve this high accuracy by averaging pion correlators over the 
location of the source operator (x,t) in Eqs. (2.1 ) and (2.5). The right panel of Fig. |2] demonstrates 
the remarkable reduction of the statistical fluctuation of C n v 4 n by this averaging. 

We determine the vector form factor Fy(q 2 ) by a constant fit to the effective value Fy(At, At';q 2 ). 
and include the correction due to the finite lattice volume estimated in one-loop ChPT [p. Although 
we do not observe any significant Q dependence of Fy(q 2 ), the spread in Fy(q 2 ) among Q = 0, —2, 
and —4 is taken as a conservative estimate of the systematic uncertainty due to the fixed topology. 

The scalar form factor normalized at a certain momentum transfer q 2 m[ can be calculated from 



F s (At,At';q 2 ) 
F s (At,At';q 2 ei ) 

R s {At,At';q 2 ) 



RsjAt^At'-q 1 ) 
R s (At,At';q 2 ei ) 

1 I 
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pMp'I fixed bl.lP'l C^(A?;p)C^(A? / ;p / ) 
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Figure 3: Left panel: effective value of normalized scalar form factor Fs(At,At';q 2 )/Fs(At,At';q 1 ef ) at 
= 0.025. Right panel: scalar form factor Fs(q 2 ) with (filled symbols) and without (open symbols) the 
contributions of the disconnected diagrams . Both data are normalized by a common value Fs(q 2 ef ) including 
the disconnected contribution. 
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At q 2 = 0, Cjisn has an additional contribution shown in Fig. [jjdue to the VEV of the scalar operator 
S. The subtraction of this contribution leads to a relatively large uncertainty in F$ (At , At' ; q 2 ) at q 2 = 
than at q 2 j^O, as seen in the left of Fig. ||| Although the Feynman-Hellmann theorem 2Fj(0) = 
dM 2 /dm ll( j provides a better determination of F s (0), it is subject to systematic uncertainties of the 
chiral extrapolation of M 2 . In our simulation setup, F$(q 2 ) has the smallest relative error at the 
smallest nonzero value of \q 2 e{ \ with |q re f| = 1. We therefore use Fs(q 2 ) normalized at this q 2 e{ in 
the following analysis. 

We determine Fs(q 2 ) /Fs(q 2 ei ) from the effective value Fs(At,At';q 2 ) /Fs(At,At';q 2 e{ ) in a sim- 
ilar way to Fy(q 2 ). The right panel of Fig. ^compares Fs(q 2 ) to that without the contribution of the 
disconnected diagrams. We observe a significant deviation between the two data, which implies 
the importance of the disconnected contributions in a precision study of F$(q ). 

4. Parametrization of q 2 dependence 

In Fig.Q, we plot the vector form factor Fy(q 2 ) and normalized scalar form factor Fs(q 2 ) /Fs(q 2 ei ) 
as a function of q 2 . We observe that Fy(q 2 ) is close to the vector meson dominance (VMD) hypoth- 
esis 1/(1 — q 2 /M 2 ) with the vector meson mass M p measured at simulated m u d. We then assume 
that the small deviation due to the higher poles or cuts can be approximated by a polynomial of q . 
The q 2 dependence of Fy(q 2 ) is therefore parametrized as 

F v(l 2 ) = i +av 1 iq 2 +ayi(q 2 ) 2 +a v ,3(q 2 ) 3 = 1 + \{r 2 ) V q 2 +c v {q 2 ) 2 + • • • (4.1) 

1 — q z /Mp o 

in order to determine the charge radius (r 2 }v and the curvature cy. This form describes our data 
well as shown in Fig. f| Results for (r 2 )y and cy do not change significantly if we remove the 
cubic term or if we add higher order terms into the parametrization form. 

Due to the lack of the knowledge about the scalar resonances at the simulated quark masses, 
we use a generic quartic form 

F s (q 2 ) = F S (0) (l + ^r 2 )^ 2 +Cs (q 2 ) 2 +a S j(q 2 ) 3 +a SA (q 2 A (4.2) 




q [GeV 2 J q~ [GeV 2 ] 

Figure 4: Vector form factor Fy(q 2 ) (left panel) and normalized scalar form factor Fs(q 2 )/Fs{q 2 e{ ) (right 
panel) as a function of q 2 . Solid lines show the fit curve and its error. We also plot the q 2 dependence of 
Fy(q 2 ) expected from the VMD model by the dashed line. 
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to parametrize the q 2 dependence of Fs(q 2 ). Our data are described by this form reasonably well as 
in Fig. f|. The result for the scalar radius (r 2 ) $ is stable against the removal of the the quartic term 
as well as inclusion of higher order terms. Such a stability is, however, not clear in the curvature 
cs due to its large statistical uncertainty. We leave a precise determination of cs for future studies, 
and only use results for (r 2 )s in the following analysis. 

ChPT can provide a more unambigu- 
ous parametrization of Fy^(q 2 )- Figure || 
shows contributions to Fy{q 2 ) from each or- 
der (q 2 ) n of a Taylor expansion of Eq. (4. 1 ). 
We observe that 0(q 6 ) and higher order con- 
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due to the use of the simple periodic bound- 4 
ary condition for quark fields. We therefore 
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pendence based on ChPT in this study. 

We note that the simulated values of 
the pion mass squared M 2 are smaller than 

(520 MeV) 2 . The 0(q 6 ) contribution to Fy(q 2 ) is small if \q 2 \ is smaller than this value. The quark 
mass dependence of our data of {r 2 )y^ is expected to be described by NNLO ChPT. 



q [GeV ] 

Figure 5: Contributions in the q 2 expansion of Fy (q 2 ) — l 
at rriud = 0.025. Thin solid, dashed and dot-dashed lines 
show 0(q 2 ), 0{q A ) and higher order contributions. The 
thick solid line is their total. 



5. Chiral extrapolation 

We first compare our lattice data of the radii (r 2 )y^ with NLO ChPT formulae [§] 



(r 2 )v 



NF 2 



(1 + 6NQ 



NF 2 



In 



M 2 



, (r 2 ) s 



NF 2 



13 

-— + 6Nf 4 



NF 2 



In 



M: 



21 



(5.1) 



where N = {An) and F is the decay constant in the chiral limit. We fix F to our estimate from 
our study of the pion decay constant [10]. The renormalization scale is set to jj. =4nF. The 



NLO fits are not quite successful as seen in Fig. |6| While our data of {r 2 )y are fitted well with 
% 2 /d.o.f. ~0.1, the value extrapolated to the physical quark mass 0.364(4) fm 2 is significantly 
smaller than experiment [12]. On the other hand, the NLO formula for {r 2 )s with the enhanced 
chiral log fails to reproduce our data and leads to large # 2 /d.o.f. ~9. 

We also note that NLO in ChPT is not sufficient to describe the quark mass dependence of 
the curvature cy. Although cy has a NLO term l/(60NF 2 M 2 ) coming from non-analytic NLO 
contributions to Fy(q 2 ), it dominates cy well below the physical pion mass and fairly near the 
chiral limit (see Fig. [7] below), where the non-analytic contributions to Fy{q 2 ) are important to be 
consistent with the existence of the 7 — > %% branch cut at q 2 >0. Since cy characterizes the 0(q 4 ) 
dependence of Fy(q 2 ), NNLO contributions are essential to describe its quark mass dependence. 

We therefore extend our analysis to NNLO ChPT. The NNLO contributions to the radii are 

P>R] 
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M 1 [GeV 2 ] 
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Figure 6: Chiral fit of (r 2 )v (left panel) and (r 2 )s (right panel) using NLO ChPT formulae. In the left panel, 
we also plot the experimental value (r 2 )v = 0.437(16) fm 2 from an analysis based on Nf = 2 ChPT [ p"T| 
(diamond) and 0.452(11) fm 2 quoted by Particle Data Group JI4 I (star). The diamond in the right panel 
represents (r 2 ),s = 0.61 (4) fm 2 obtained from an indirect determination through %% scattering 
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The NNLO expression of cy is 
1 1 



Cy 
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8429 7V„. AT Ar2 r 



6(WF 2 M 2 A^f 4 ^720 25920 3 M 6 
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(5.2) 



(5.3) 



(5.4) 



The analytic terms with T r sy S \ r c ~, represent contributions from 0(p 6 ) chiral Lagrangian. The linear 
combination l\— 1^/2 appearing commonly in (r 2 )y and cy is denoted by l[ 2 for simplicity. 




7 



Pion form factors from lattice QCD with exact chiral symmetry 




We first carry out a simultaneous fit to (r 2 }y and cy in terms of M 2 /(NF 2 ). This fit has only 
four free parameters Z£, l[ 2 > r vr anc ^ r v c !01> an d these can be determined with reasonable accuracy 
without introducing phenomenological inputs. As seen in Fig. 0, our data are well described by the 
NNLO formulae with # 2 /d.o.f. ~0.7. The extrapolated values of (r )y and cy are consistent with 



recent phenomenological determinations [11, |15J, [16 ]. 

The inclusion of (r 2 )s into the simultaneous fit introduces additional four free parameters 
V 2 , I3, l\ and r r Sr , and we need to fix some of them to obtain a stable fit. In this study, we use a 
phenomenological estimate 1% =4.31(1 1) Jl3| ] and a lattice estimate % = 3.38(56) from our analysis 
of the pion mass JTo[ ] 1 , since they are determined with a reasonable accuracy and appear only in 
the NNLO terms. As plotted in Fig. |8[ this fit describes our data of (r 2 )s reasonably well with 
^ 2 /d.o.f. ~0.7. The extrapolations of (r 2 )v and cy are consistent with those in Fig. [7]. 

At physical quark mass, we obtain 

(r 2 ) v = 0.409(23)(37) fm 2 , (r 2 ) s = 0.617(79)(66) fm 2 , c v = 3.22(17)(36) GeV 4 , (5.5) 



where the first and second errors are statistical and systematic, respectively. The latter includes 
uncertainties due to the choice of the input to fix the lattice scale and the inputs for the LECs (F, 
h and h) as well as uncertainties due to the chiral extrapolation and lattice discretization. These 
results for (r 2 )y,s and cy are consistent with phenomenological analyses. 
The results for the relevant LECs are 



f 6 = 11.9(0.7)(1.0), U= 4.09(50) (52), h -l 2 = -2.9(0.9) (1.3), (5.6) 
r^ = -1.0(1.0)(2.5)xl0- 5 , r r Vc 



1.0(1.0)(2.5) x 10~ 5 , r^ c = 4.00(17) (64) x 10~ 5 , r' Sr = 1.74(36) (78) x 10~ 4 . (5.7) 



Our estimate of l(, is slightly smaller than those from ChPT analyses: l(, = 16.0(0.9) from Fy [11] 
and 15.2(0.4) from T and 71 decays [O]. This is partly because our estimate of F JlOt ] is slightly 
smaller than phenomenological estimates. We note that I4 is consistent with our determination 
I4 =4.12(56) fromF^ [ ^0| ] and a phenomenological estimate 4.39(22) [13]. 



'The fl- independent convention h, is defined by l\ = ji(Ji + \v\M\j 'jl 2 ])/2N with y\ = 1/3, Ji = 2/3, f} = — 1/2, 
74 = 2 and y6 = -l/3. 
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6. Conclusions 

In this article, we present our calculation of pion form factors in two-flavor lattice QCD with 
exact chiral symmetry, which enables us to unambiguously compare our lattice data with two-loop 
ChPT. By employing the all-to-all quark propagators, Fs(q 2 ) is calculated including contributions 
from the disconnected diagrams for the first time. We observe that two-loop contributions are 
important to describe the quark mass dependence of (r)v,s and cy at our region of the pion mass 
M K > 290 MeV. Our chiral extrapolation of (r 2 )^ and cy are consistent with phenomenological 
analyses. We also confirm that Fs(q 2 ) and F % lead to consistent results for I4. 

For a more precise comparison with experiment, we need to extend this study to three-flavor 
QCD. Such simulations are currently underway. Another important subject is a better control of 
the parametrization of the q 2 dependence of Fy } s(q 2 )- To this end, the use of the twisted boundary 
condition [18] to simulate small values of |<7 2 |, dispersive analysis of the q 2 dependence [[TJ[ |TJ], 



and model-independent determination of the scalar resonance mass [J19J are interesting possibilities 
for our future studies. 
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